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ABSTRACT 

Using a complete sample of ~300 star-forming galaxies within 11 Mpc of the 
Milky Way, we evaluate the consistency between star formation rates (SFRs) 
inferred from the far ultraviolet (FUV) non-ionizing continuum and Ha nebular 
emission, assuming standard conversion recipes in which the SFR scales linearly 
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with luminosity at a given wavelength. Our analysis probes SFRs over 5 orders of 
magnitude, down to ultra-low activities on the order of ~10~ 4 M yr _1 . The data 
are drawn from the 11 Mpc Ha and Ultraviolet Galaxy Survey (11HUGS), which 
has obtained Ha fluxes from ground-based narrowband imaging, and UV fluxes 
from imaging with GALEX. For normal spiral galaxies (SFR~1 M Q yr _1 ), our 
results are consistent with previous work which has shown that FUV SFRs tend 
to be lower than Ha SFRs before accounting for internal dust attenuation, but 
that there is relative consistency between the two tracers after proper corrections 
are applied. However, a puzzle is encountered at the faint end of the luminosity 
function. As lower luminosity dwarf galaxies, roughly less active than the Small 
Magellanic Cloud, are examined, Ha tends to increasingly under-predict the total 
SFR relative to the FUV. The trend is evident prior to corrections for dust 
attenuation, which affects the FUV more than the nebular Halpha emission, so 
this general conclusion is robust to the effects of dust. Although past studies 
have suggested similar trends, this is the first time this effect is probed with a 
statistical sample for galaxies with SFR<0.1 M Q yr" 1 . By SFR~0.003 M yr" 1 , 
the average Ha-to-FUV flux ratio is lower than expected by a factor of two, and 
at the lowest SFRs probed, the ratio exhibits an order of magnitude discrepancy 
for the handful of galaxies that remain in the sample. A range of standard 
explanations does not appear to be able to account for the magnitude of the 
systematic. Some recent work has argued for an IMF which is deficient in high 
mass stars in dwarf and low surface brightness galaxies, and we also consider this 
scenario. Under the assumption that the FUV traces the SFR in dwarf galaxies 
more robustly, the prescription relating Ha luminosity to SFR is re-calibrated 
for use in the low SFR regime when FUV data are not available. 

Subject headings: galaxies: evolution — galaxies: dwarf — ultraviolet: galaxies 
— galaxies: photometry — stars: formation — surveys 

1. Introduction 

Ha nebular line emission and the ultraviolet (UV) non-ionizing continuum flux are 
two fundamental, widely-used star formation rate (SFR) indicators. Ha nebular emission 
arises from the recombination of gas ionized by the most massive O- and early-type B-stars 
(M* > 17M ). It therefore traces star formation over the lifetimes of these stars, which is 
on the order of a few million years. In contrast, the UV flux primarily originates from the 
photospheres of a fuller mass spectrum of O- through later-type B-stars (M* > 3M ), and 
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thus measures star formation averaged over a longer ~ 10 8 yr timescale. 

In principle, the two tracers should yield consistent SFRs. Discrepancies that are other- 
wise found can be valuable in revealing problems in the calibrations of the indicators that are 
used and/or the assumptions underlying them. A number of studies have directly compared 
the integrated Ha and UV SFRs of galaxies for which both diagnostics have been measured 
(e.g., Buat et al. 1987; Buat 1992; Sullivan et al. 2000; Bell & Kennicutt 2001; Buat et 
al. 2002; Iglesias- Paramo et al. 2004; Salim et al. 2007), with most assuming standard 
conversion recipes where the ratios of the UV and Ha luminosities to the SFR are constants 
(e.g. Kennicutt 1998; hereafter K98). For normal spiral galaxies (SFR>1 M yr -1 ), there 
is broad agreement among such studies that although the UV SFR tends to be lower than 
the Ha SFR before accounting for internal dust attenuation, there is good consistency be- 
tween the two after proper corrections have been applied. This implies that the UV stellar 
continuum is more affected by attenuation than the Ha nebular emission, and the intrinsic 
dust-corrected Ha-to-UV flux ratio is constant on average. As lower mass dwarf galaxies are 
probed however, there have been some indications that the observed Ha-to-UV flux ratio 
may systematically decrease, and SFRs inferred from Ha will tend to be lower relative to the 
UV (e.g., Sullivan et al. 2000; Bell & Kennicutt 2001). Such a trend could provide evidence, 
for example, that the recent star formation histories of low-mass galaxies are predominantly 
bursty, or otherwise exhibit abrupt changes. Following a discontinuity in star formation, 
a dearth of ionizing stars develops as those stars expire relative to the longer-lived, lower 
mass stars that also significantly contribute to the UV emission, and this may cause the 
Ha-to-UV flux ratio to be lower than expected (e.g., Sullivan et al. 2004; Iglesias-Paramo et 
al. 2004). Or it may perhaps suggest that the stellar initial mass function is not universal 
and becomes more deficient in the highest mass ionizing stars within low-mass, low-surface 
brightness galaxies (e.g., Meurer et al. 2009; Pflamm-Altenburg et al. 2009). However, 
wholesale deviations from the expected ratio as a function of luminosity have only begun 
to become apparent when integrated SFRs less than about 0.1 M yr -1 have been probed 
(i.e., a few times lower than the SFR of the Large Magellanic Cloud), and only handfuls 
of galaxies with such low activities were included in the previous studies. Therefore, the 
observational trend first needs to be more robustly established, preferably with an unbiased, 
statistical sample of dwarf galaxies. 

In order to further study the relationship between the UV and Ha emission, we draw 
upon GALEX UV observations for a complete sub-set of star- forming galaxies from the 11 
Mpc Ha narrowband imaging survey of Kennicutt et al. (2008). Dwarf galaxies dominate this 
sample by number, and over 80% have Ha-based SFRs<0.1 M yr -1 . This dataset is thus 
ideal for probing the relationship between the UV and Ha emission in the low SFR regime. 
Such a sample is also particularly well-suited for study in the UV since the vast majority of 



-4- 



the galaxies have low-metallicities and low-dust contents which minimizes difficulties with 
otherwise vexing attenuation corrections. 

The precursor Ha survey and GALEX follow-up program together has been referred 
to as 11HUGS, the 11 Mpc Ha and Ultraviolet Galaxy Survey. The data from 11HUGS 
are part of a broader campaign designed to provide a census of star formation and dust in 
the Local Volume. 11HUGS is being further augmented by Spitzer IRAC mid-infrared and 
MIPS far-infrared observations through the composite Local Volume Legac)Q (LVL) program 
(Lee et al. 2008; Dale et al. 2009), which is delivering the overall multi-wavelength dataset 
to the community!^ The ensemble of observations enables studies ranging from detailed 
spatially-resolved analyses probing the properties of individual HII regions and the radial 
profiles of disks (e.g., Munoz-Mateos et al. 2009), to a broader statistical examination of 
the demographics of the local star- forming galaxy population as a whole (e.g., Lee et al. 
2007, Dale et al. 2009). Other specific studies currently being pursued by the 11 HUGS and 
LVL teams include the calibration of monochromatic far-IR SFR indicators (Calzetti et al. 
2009) and photometric measurement indicies for PAH emission (Marble et al. 2009), and 
an examination of the prevalence and duty cycle of starbursts in dwarf galaxies (Lee et al. 
2009a). 

An outline of this paper is as follows. In Section 2, we describe the 11HUGS Ha nar- 
rowband and GALEX UV datasets which are used in this analysis. In Section 3, we compute 
SFRs and evaluate the consistency between the two tracers over 5 orders of magnitude down 
to ultra-low activities of 10~ 4 M yr -1 . We adopt the SFR conversion prescriptions of K98, 
so this is equivalent to examining whether the Ha-to-UV flux ratio exhibits systematic de- 
viations from a constant value. Our primary observational result is that the Ha-to-FUV 
flux ratio decreases with decreasing B-band, Ha and FUV integrated luminosities. This 
systematic is evident prior to any corrections for internal dust attenuation, which affects the 
FUV more than the Ha nebular emission. Our analysis provides a robust confirmation of 
trends suggested by prior studies, with extension to lower SFRs than previously probed. In 
Section 4, we investigate possible causes of this trend. The examination of the assumptions 
underlying the K98 calibration provides a convenient framework for our analysis, and we 
check the validity of these assumptions with respect to the local dwarf galaxy population. 
We discuss the effects of uncertainties in the stellar evolution tracks and model atmospheres, 
non-solar metallicities, non-constant SFHs, possible leakage of ionizing photons/departures 



x http: / /www. ast.cam.ac.uk/research/lvls/ 

2 Public data releases have begun through the NASA/IPAC Infrared Science Archive (IRSA; 
http://ssc.spitzer.caltech.edu/legacy/lvlhistory.html). 
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from Case B recombination, dust attenuation and stochasticity in the formation of high-mass 
stars. However, none of these drivers, at least when considered individually, appear to be 
able to account for the observed trend. We therefore also consider systematic variations in 
the IMF where dwarf galaxies with lower activities are deficient in the most massive ionizing 
stars. At the end of Section 4, we re-evaluate Ha as a SFR indicator for dwarf galaxies, and 
provide an alternate conversion prescription for use when FUV luminosities are not available. 
We conclude in Section 5 by summarizing our results, and outline remaining uncertainties 
in the interpretation and possible directions for future work. 

2. Data 

To investigate the relationship between the FUV and Ha emission in the low luminosity 
regime, we draw upon data collected by the 11HUGS GALEX Legacy program. The 11HUGS 
sample is dominated by dwarf galaxies and is thus ideal for studying the nature of galaxies 
with low SFRs. Integrated Ha and UV flux catalogs are provided by Kennicutt et al. (2008; 
hereafter Paper I) and Lee et al. (2009b; hereafter Paper II), respectively. Details on the 
sample selection, observations, photometry and general properties of the sample are given 
in those papers. Here, we provide a brief summary of the dataset. 

Our overall Local Volume sample of 436 objects is compiled from existing galaxy catalogs 
(as described in Paper I), and is divided into two components. The primary component aims 
to be as complete as possible in its inclusion of known nearby star-forming galaxies within 
given limits, and consists of spirals and irregulars (T>0) within 11 Mpc that avoid the 
Galactic plane (\b\ > 20°) and are brighter than B — 15 mag. These bounds represent the 
ranges within which the original surveys that have provided the bulk of our knowledge on the 
Local Volume galaxy population have been shown to be relatively complete (e.g., Tully 1988, 
deVaucouleurs et al. 1991), while still spanning a large enough volume to probe a diverse 
cross-section of star formation properties. A secondary component consists of galaxies that 
are still within 11 Mpc, but fall outside one of the limits on brightness, Galactic latitude, or 
morphological type, and have available Ha flux measurements (i.e., targets that were either 
observed by our group as telescope time allowed, or had Ha fluxes published in the literature). 
Subsequent statistical tests, as functions of the compiled B-band apparent magnitudes and 
21-cm (atomic hydrogen, HI) flux, show that the overall sample of star-forming galaxies is 
complete to ~15.5 mag and 6 Jy km s _1 , respectively (Lee et al. 2009a). These limits 
correspond to M B < -15 and M HI > 2 x 10 8 M for |6| > 20° at the edge of the 11 Mpc 
volume. 

Through a combination of new narrowband Ha + [Nil] and i?-band imaging, and data 
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compiled from the literature, Paper I provides integrated Ha fluxes for over 90% of the 
total sample. The new narrowband data were obtained at 1-2 m class telescopes in both 
hemispheres, and reached relatively deep point source flux and surface brightness limits of 
~2 x 10~ 16 ergs cm~ 2 s _1 and ~4 x 10 -18 ergs cm -2 s _1 arcsec -2 , respectively. 

The GALEX UV follow-up imaging primarily targeted the |6| > 30°, B < 15.5 subset 
of galaxies. The more restrictive latitude limit was imposed to avoid excessive Galactic 
extinction and fields with bright foreground stars and/or high background levels for which 
UV imaging would be prohibited due to the instrument's brightness safety limits. The 
GALEX satellite utilizes a 50 cm aperture telescope with a dichroic beam splitter that 
enables simultaneous observations in the FUV (A e //=1516A, FWHM=269A) and the NUV 
(A e// =2267A, FWHM=616A). Martin et al. (2005), Morrissey et al. (2005), Morrissey et al. 
(2007) provide full details on the satellite, telescope, instrument and calibration and data 
processing pipeline. Deep, single orbit (~1500 sec) imaging was obtained for each galaxy, 
following the strategy of the GALEX Nearby Galaxy Survey (Gil de Paz et al. 2007). 
GALEX observations for a significant fraction of the remaining galaxies beyond these limits 
have been taken by other programs and are publicly available through MAST§] Paper II 
provides photometry (following the methods of Gil de Paz et al. 2007) based upon both 
the archival data and 11HUGS imaging for over 85% of the overall Paper I sample, and the 
currently available measurements are incorporated into our analysis. Table 1 lists the ~300 
galaxies included in this study. 



3. Comparison of Ha and FUV SFRs 

To begin, we compute the FUV and Ha SFRs, where we correct for Milky Way fore- 
ground reddening but not for attenuation internal to the galaxies. The standard conversion 
prescriptions of K98, where the ratio of the luminosities to the SFR are constants, are as- 
sumed. We first describe the calculation of the SFRs, and then examine the consistency 
between the two tracers. 

Corrections for foreground reddening are calculated by using E(B — V) values based on 
the maps of Schlegel et al. (1998) and the Cardelli et al. (1989) extinction law with Ry = 3.1. 
The resulting relationships between the color excess and extinction are A FUV = 7.9E(B — V) 
and A Ha = 2.5E(B-V). 



3 MAST is the Multimission Archive at the Space Telescope Science Institute 
( http: / / archive.stsci.edu/index.htnu] ) . 
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We compute Ha SFRs using the emission-line fluxes given in Table 3 of Paper I. The 
narrowband filters used for the survey are wide enough to capture the [NII]AA6548,83 lines 
which flank Ha, so the [Nil] component of the observed flux needs to be removed before SFRs 
can be calculated. Measurements of [Nil] /Ha from spectroscopy are available for about a 
quarter of the sample and are used to correct the fluxes in those cases. A compilation of 
measurements from the literature is provided in Table 3 of Paper I. Otherwise, the ratio is 
estimated from the average relationship between [NII]A6583/Ha and M B , a consequence of 
the luminosity-metallicity relationship for galaxies. The adopted scaling relationship, derived 
using the integrated spectroscopic dataset of Moustakas & Kennicutt (2006), is: 

log([NII]A6583/Ha) = (-0.173 ± 0.007) M B - (3.903 ± -0.137) if M B > -20.3 

[NII]A6583/Ha = 0.54 if M B < -20.3 ^ ' 

The total [NII]AA6548,83/Ha ratio is computed assuming a 3-to-l ratio between [NII]A6583 
and [NII]A6548 (Osterbrock 1989). Further details can be found in Appendix B of Paper 
I. As noted there, although estimates based on this correlation are only good to a factor of 
~2, the majority (>70%) of galaxies in our dwarf-dominated sample have M B > —17, and 
hence have [NII]/Ha ratios <0.1. The corresponding uncertainties in the Ha flux due to 
[Nil] corrections are thus lower than <10%. 

FUV SFRs are based on the FUV integrated photometry in Paper II, Table 2. /„ 
follows directly from the tabulated AB magnitudes given there since /^[ergs s" 1 cm -2 Hz -1 ] 
= io~°- 4(mAS + 48 - 6 ). Calculations based on the NUV data will generally be less reliable for 
tracing the recent star formation since the flux at these redder NUV wavelengths will have 
a greater contribution from stars with lifetimes >100 Myrs. 

FUV and Ha luminosities then follow from the distances in Paper I, Table 1, and the 
SFRs are computed as in K98: 

SFR[M yr^] = 7.9 x 10~ 42 L(Ha) [ergs s" 1 ] (2) 
SFRfMoyr" 1 ] = 1.4 x 10~ 28 L„(UV) [ergs s^Hz" 1 ] (3) 

The K98 conversion factors are calculated using a Salpeter IMF with mass limits of 0.1 
and 100 M , and stellar population models with solar metal abundance. The Ha and UV 
calibrations also assume that the star formation history is constant for at least the past ~10 
Myrs and ~100 Myrs, respectively. Table 1 provides the computed SFRs along with some 
general properties of the galaxies in the sample. The assumptions underlying the conversions 
and their validity in the context of the our dwarf galaxy dominated Local Volume sample 
are discussed further in the next section. 
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In Figure [H we plot these non-dust-corrected FUV SFRs against the Ha SFRs. The solid 
line represents a one-to-one correlation between the SFRs. Axes indicating the corresponding 
luminosities are also shown. Different colors and symbols are used to distinguish between 
morphological types as specified in the figure. The error bars reflect uncertainties in the 
flux measurements, as described in Papers I and II. A wide range of integrated SFRs are 
covered — the sample probes normal spiral disks with SFRs of about unity to dwarf irregular 
galaxies with ultra-low SFRs of ~10~ 4 M yr _1 . 

In general, there is coarse agreement between the FUV and Ha SFRsE The SFRs 
agree to within about a factor of 2 for the majority of galaxies with SFR>0.01 M yr _1 , 
even before internal extinction corrections are applied. This agreement is reasonable given 
the relative absence of highly-obscured, luminous galaxies in the small volume probed by 
our sample. However, there is a clear systematic present: at low SFRs (log SFR~-1.5; 
SFR<0.03M Q yr _1 ) Ha begins to underpredict the SFR relative to the FUV luminosity. A 
linear least squares bisector fit to the data demonstrates this tilt (dashed linejE 

The systematic offset is better illustrated in the top panel of Figure [2] which plots the 
ratio of the SFRs as a function of the Ha SFR. The covariance between the axes acts to 
emphasize the deviation from unity. The right and top axes indicate the corresponding 
flux ratios and luminosities, respectively. Here, a linear least squares fit to the data is 
shown (dashed line) and has a scatter of 0.30 dex. Binned averages with la ranges are also 
overplotted (blue symbols), and the values are given in Table [2J For L(Ha)> 10 40 ergs s _1 
(SFR> 0.1M Q yr" 1 ) there is good consistency between the two tracers, and the average 
deviation from unity is 30%. For L(Ha) < 4 x 10 38 ergs s" 1 (SFR< O.OO3M yr" 1 ), Ha 
yields SFRs that are lower than the FUV by average factors ranging from two to more than 
ten. An analogous plot as a function of the B-band absolute magnitude is shown in the 
bottom panel of Figure [2J Here, the tilt in the correlation is less steep, and the scatter is 
larger (0.35 dex). For galaxies brighter than —19 mag the Ha yields an SFR that is ~50% 
higher. There is good consistency for —16 > Mb > —19, and average offsets of a factor of 2 



4 There is one significant outlier at low luminosities where the Ha SFR appears to exceed the FUV SFR 
by more than a factor of 10. This system, KDG61, is in the M81 group, and Croxall et al. (2009) speculate 
that the system is a chance superposition of a dwarf spheroidal galaxy and an HII region from the M81 tidal 
stream. Such a scenario can plausibly result in a high Ha-to-UV ratio and the extremely red color (FUV- 
NUV>2, which is not the result of extinction since Ha/H/3 is close to the expected Case B recombination 
value) that are observed. This data point is excluded from the least squares fits that follow. 

5 The filled black symbols represent galaxies where the Ha flux may be underestimated because the 
narrowband imaging did not wholly enclose the galaxy. These points are excluded from the fit, and omitted 
from Table 1. In the next figure, these same galaxies are plotted with the circled points, and are also excluded 
from the fit and reported statistics. 
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and higher begin to appear at > — 14 mag. 

4. Discussion: Understanding the Systematic Decline in L(Ha)/L(FUV) 

Using a complete, statistical sample of star-forming galaxies within the Local Volume, 
we find that the observed Ha-to-FUV flux ratio systematically decreases with declining SFR 
and galaxy luminosity. When standard conversion recipes of the form SFR/L=constant are 
applied (e.g., K98), the Ha luminosity will thus underestimate the SFR relative to the FUV 
luminosity in dwarf galaxies. This analysis confirms earlier indications of a such a trend 
for L(Ha)< 10 40 or SFR<0.1 M yr _1 (see introduction for references). Our more robust 
sampling of low-mass star-forming galaxies represents a significant improvement over past 
studies which were only able to probe to regimes of —2 < log SFR(Ha)< —1 with ~20 
dwarf galaxies (c.f. Sullivan et al. 2000, Figure 14; Bell & Kennicutt 2001, Figure 2). The 
recent study of Meurer et al (2009), which has examined correlations in the Ha-to-FUV flux 
ratio with Ha and i?-band surface brightness, also only includes ~25 galaxies with —3 < 
log SFR(Ha)< —1 (M. Seibert, private communication). In contrast, the current analysis is 
based on over 200 dwarf galaxies with SFR<0.1 M yr _1 and extends to ultra-low activities 
of SFR~ 10" 4 M yr" 1 . We find that Ha SFRs of 3 xl0~ 3 M Q yr" 1 are lower than those 
inferred from the FUV by a factor of ~2 and the discrepancy increases to a factor of >10 
by SFR(Ha)=10~ 4 M yr -1 , on average. 

The expectation that the two tracers should yield consistent SFRs relies on the robust- 
ness of the SFR calibrations that are used. In the K98 formulation, consistent SFRs require 
that L(Ha)/L(FUV)=constant, since SFR/L(Ha) and SFR/L(FUV) are also both simple 
constants. Therefore, a convenient framework for understanding the observed deviations 
from the expected ratio involves checking the validity of the assumptions underlying the cal- 
ibrations with respect to the particular galaxies to which they are being applied. To briefly 
summarize how the constants are derived, stellar population synthesis models are used to 
compute Ha and UV luminosities, given an invariant stellar initial mass function (IMF). The 
expected flux is calculated by integrating a constant SFR over a time period that is long 
enough to allow the birth and death of the stars responsible for the emission to equilibrate. 
Solar metallicity models are normally used, and the IMF is well-populated at all masses^ 



6 For example, the Starburst99 models of Leitherer et al. (1999) use 10 6 M© single-age stellar populations 
to construct systems with more complex star formation histories, and such building blocks contain ^4500 
O-stars at zero-age for a Salpeter IMF with mass limits of 1 and 100 Mq yr _1 . For more commonly used 
mass limits of 0.1 and 100 M© yr -1 , the corresponding number of O-stars would be lower, but still well 
sampled, at ^2500. 



-10- 



The nebular emission luminosity is computed assuming Case B recombination to convert 
the production rate of ionizing photons to an Ha luminosity. Finally, the model-predicted 
luminosities are typically intrinsic values that have not yet been attenuated by dust. 

All of these ingredients have been examined in previous studies to investigate possible 
trends in the L(Ha)/L(FUV) ratio, but primarily in the context of more massive star-forming 
systems than typical of our Local Volume sample. The effects of uncertainties in the stellar 
evolution tracks and model atmospheres, non-solar metallicities, non-constant SFHs, leakage 
of ionizing photons, departures from Case B recombination, dust attenuation, stochasticity in 
the formation of high-mass stars, and finally, variations in the IMF have all been considered 
(e.g., Buat et al. 1987; Buat 1992; Sullivan et al. 2000; Bell & Kennicutt 2001; Buat et 
al. 2002; Iglesias-Paramo et al. 2004; Sullivan et al. 2004; Salim et al. 2007) and have 
been recently revisited by Meurer et al. (2009). More mundane measurement systematics 
such as photometric aperture mismatches between the FUV and Ha data have also been 
considered^ (e.g., Sullivan et al. 2000). We re-evaluate these issues in the context of our 
dwarf-dominated sample, and tackle the ones which are more straightforward to examine 
first. 



4.1. Internal Dust Attenuation 

In general, attenuation by dust is the first issue to address when comparing UV/optical 
observed luminosities with intrinsic model values. For normal spiral galaxies (SFR>1 M Q yr _1 ) 
there is broad agreement among previous work that although the UV SFR tends to be lower 
than the Ha SFR before accounting for internal dust attenuation, there is good consistency 
between the two after proper corrections have been applied. Empirically, these results sug- 
gest that the UV stellar continuum is more affected by attenuation than the Ha nebular 
emission^]. If the same holds true for lower mass dwarf galaxies, dust corrections would drive 



7 The Ha fluxes given in Paper I and the FUV magnitudes given in Paper II that are used in the SFR 
comparison analysis represent integrated measurements which are not necessarily taken through the same 
apertures. However, we have checked that the impact of this on our analysis and conclusions is negligible. 
This is because the GALEX and Ha apertures cover the same general area for most galaxies. There are 
some cases, however, where the Ha aperture is smaller because (i) it is limited by the FOV of the optical 
detector or (ii) the FUV emission is more extended. In such cases, GALEX fluxes were re-measured using 
the apertures used for the H-alpha photometry and there are no significant differences in the correlations 
when these fluxes are used. 

8 Based on extinction curves alone, it may seem obvious that the UV should be more attenuated than 
redder emission at 6563A. However, differential reddening between gas and stars must also be considered. 
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the Ha-to-UV flux ratio even lower, which would at first appearance seem to exacerbate 
rather than help reconcile the discrepancy. To check this for our sample, we estimate the 
attenuation in the Ha flux with the Balmer decrement, and in the FUV flux using the to- 
tal infrared-to-UV (TIR/FUV) flux ratio. The two methods are independent, the former 
being rooted in recombination physics and the latter in stellar synthesis and dust models, 
although an extinction curve must be assumed in each case. The resulting dust corrections 
are provided in Table 1. 

4-1.1. via the Balmer Decrement 

The validity of using the Balmer decrement to infer the average amount of nebular 
extinction in galaxies is supported by studies that have shown that the scatter between 
infrared-based SFRs and Ho; SFRs is decreased when the Ha luminosity is corrected for 
extinction in this manner (e.g., Dopita et al. 2002; Kewley et al. 2002; Moustakas et al. 
2006). Kennicutt et al. (2009) reaches similar conclusions from the analysis of about 400 
nearby galaxies. As in Lee (2006), we assume an intrinsic Case B recombination ratio of 
2.86 for Ha/H/5, and the Cardelli et al. (1989) Milky Way extinction curve with R v = 3.1, 
to express An a in terms of the observed Ha/H/5 ratio as: 

A Ha = 5.91 log^ - 2.70. (4) 

JH/3 

For ~20% of the sample, spectroscopic measurements of H.a/H/3 from the literature are 
available. For those without robust spectroscopic measurements, we follow a strategy similar 
to that used for estimating [Nil] /Ha above. Using the integrated spectroscopy of Moustakas 
& Kennicutt (2006) we derive an empirical scaling relation between A^ a and Mb which is 
shown in Figure [3J The data exhibit the well-known trend that the more luminous galaxies 
are more heavily obscured. We fit a piecewise function to the data, described by a constant 
value at the lowest luminosities and a second order polynomial for Mb < —15.0. Objects 
that have An a < 0, a signal-to-noise less than 10 in the H/3 line, an H/3 equivalent width 
< 5 A, or show evidence of AGN emission are excluded from the fit. The adopted scaling 
relation is: 



As discussed in Calzetti et al. (1994), the nebular emission (from HII regions which are the sites of the 
most recent star formation and are enshrouded by dust) suffers about twice the reddening experienced by 
the stellar continuum, and this leads to an Ha attenuation that is more comparable to that of the UV 
continuum. 
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f 0.10 if M B > -14.5 

^ Q \ 1.971 + 0.323Mb + 0.0134M| ifM B <-14.5 U 

Of course, there is a great deal of scatter from this average relationship which becomes 
more severe with increasing luminosity. However, although there is a 40% scatter for galaxies 
with Mb < —18, the scatter is much lower for the relatively transparent dwarf galaxies which 
dominate our Local Volume sample: for —14.7 > Mg > —18 it is 20% and for the lowest 
luminosity galaxies, where we have assumed a constant average correction, it drops further 
to 10%. 



4-1-2. A FU y 

The total infrared-to-UV flux ratio (TIR/FUV) is an indicator of the UV attenuation 
in a galaxy as it represents the total stellar light which has been absorbed and re-radiated 
by dust, relative to the surviving stellar UV light which is directly observed. Using stellar 
population models with various dust geometries and extinction curves, it has been shown 
that for systems with recent star-formation, such as the spiral and irregular galaxies in our 
sample, the mapping between A(FUV) and TIR/FUV is relatively tight (~ 20%); that is, 
the contribution of less massive/older stellar populations to the heating of the dust is either 
negligible or is in direct proportion to that of the O&B stars which dominate the UV emission 
(e.g., Buat et al. 2005). The mapping as given by Buat et al. is: 



A(FUV) = -0.0333x 3 + 0.3522x 2 + 1.1960x 3 + 0.4967 (6) 

where x=log(-|^). MIPS 24, 70 and 160/im integrated fluxes from the Spitzer LVL program 
are available for the majority of galaxies used in this analysis (Dale et al. 2009), and we 
compute TIR fluxes for those galaxies that have robust detections in all three bands (50% 
of our sample) using the calibration of Dale & Helou (2002): 



TIR = 1.5591/^(24) + 0.7686^(70) + 1.347^(24). (7) 

FUV is computed as vj v at 1520A. 

For those galaxies for which both a spectroscopic measurement of Ha/H/5 and MIPS FIR 
data are available, we compare the Ha and FUV attenuations in bottom panel of Figure H] 
(solid points). There is a good correlation with a slope of 1.8 (solid line). Interestingly, 
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A(FUV)/A(Hot)=1.8 is what is expected for the Calzetti obscuration curve and differential 
extinction law (Calzetti 2001), which is based on UV luminous starbursts. This agreement 
provides some assurance that the attenuation corrections we have derived are reasonable and 
generally consistent. In Figure H] we also plot the A(Ha) estimated from Mb using equation 
[5] (open circles). As expected, there is more scatter between A(FUV) and the estimated 
A(Ha), but the values are still correlated with a slope of ~1.8. 

When FIR data are not available, it is common to resort to methods which estimate 
A(FUV) directly from some measure of the UV spectral slope (e.g., a UV color). Previously, 
this method appeared to be reasonable since it is well-known that the UV spectral slope j3 
and TIR/FUV are tightly correlated (and thus the UV slope and A(FUV) should also be 
well-correlated), at least for samples of local starbursts (Calzetti et al. 1994; Meurer et al. 
1999). It has become clear however that the "IRX-/3" correlation for starbursts does not 
hold for normal star-forming galaxies, which have lower values of TIR/FIR at a given UV 
color, and more scatter between the quantities. Calibrations and scaling relationships which 
better describe the correlation for such systems have been published (e.g., Kong et al. 2004; 
Seibert et al. 2005; Cortese et al. 2006; Boisser et al. 2007; Gil de Paz et al. 2007). However, 
as we show in Dale et al. (2009) and Lee et al. (2009b), the majority of galaxies in our local 
volume sample are low-luminosity dwarf irregulars which are blue (0<FUV-NUV<0.5), have 
low TIR/FUV (< 2), and thus lie in the area of the IRX-/3 diagram where the UV color 
is not (strongly) correlated with TIR/FUV and attenuation. Thus, such IRX-/3 relations 
cannot be used to infer TIR/FUV and A(FUV) for our sample overall. For example, a 
relation by Cortese et al. (2006), which is based on normal but more luminous spirals and 
provides a poor fit for galaxies with TIR/FUV< 2 results in A(FUV) estimates that are too 
high when compared to the attenuations based on the Balmer decrement, as illustrated in 
Figure HI We thus take another approach for estimating A(FUV) for galaxies with no FIR 
data or only upper-limit detections. Based on Figure HJ we simply scale the computed A (Ha) 
values (whether calculated from the Balmer decrement or estimated from Mb) by a factor of 
1.8. This value is also adopted when Equations [6] and [7] produce a negative correction (i.e., 
for TIR/FUV<0.3). In all, A(FUV) is estimated from A(Ha) for ~60% of the galaxies in 
Table 1. Based on comparison with TIR/FUV based attenuations (Figure H]), there is a lcr 
uncertainty of ~17% when estimating A(FUV) from Balmer decrement based A(Ha) values. 
When A(Ha) is itself estimated (equation 5), the average uncertainties increase (26%), and 
are larger for more luminous objects {Mb < —18; ~40%) than for the dwarf galaxies that 
dominate our local volume sample {Mb > —18; ~20%). 
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4-1.3. Comparison of Dust Corrected SFRs 

After accounting for the effects of dust using the best available attenuation value for 
each galaxy, the Ha-to-FUV flux and SFR ratios are replotted in Figure [51 Data points 
where more robust corrections are applied (i.e., based on Balmer decrements and/or TIR- 
to-FUV ratios) are distinguished in red. There are no significant differences between the 
trends described by galaxies where scaling relationships are used to estimate attenuations 
and those which have more robust corrections. Binned averages with la scatters are again 
given in Table [2j For all galaxies, the correction increases the FUV flux relative to the Ha 
flux (FUV is more attenuated as expected). However, since higher luminosity galaxies tend 
to suffer from more attenuation than those at lower luminosity (e.g., Figure [3]), the Ha-to- 
FUV flux ratio is depressed by a greater factor at the high luminosity end. At L(Ha)~ 10 41 
ergs s _1 (SFR ~ 1M Q yr _1 ), the ratio decreases by ~ 0.2 dex, and as a result falls below 
the K98 value (solid line). Galaxies with L(Ha)< 10 38 ergs s" 1 (SFR(Ha) < 1CT 3 M yr" 1 ) 
are minimally affected. The main consequence is that the slope above SFR(Ha) ~ 3 x 10 -2 
M Q yr _1 is flattened (i.e., the dust-corrected ratio in this regime is constant on average). 
Adopting this lower value as the fiducial expected ratio (instead of the K98 value) would 
have the effect of mitigating the relative discrepancy at lower luminosities by 0.1 dex. Even 
in this situation however, factor of 2 offsets in the Ha-to-FUV ratio at SFR(Ha) ~ 2 x 10~ 3 
M yr -1 , which increase to factors of > 10 at SFR(Ha)=10~ 4 M yr" 1 , would still remain. 

4.2. Stellar Model Uncertainties 

Differences in stellar evolution and atmosphere models used to calibrate Ha and FUV 
luminosities as SFR indicators give rise to differences in the respective SFR conversion factors 
and hence to the expected Ha-to-FUV ratio. While this will not produce systematic trends 
in the ratio as a function of the luminosity for a given metallicity and IMF, it does define the 
fiducial from which deviations are measured. In Figure a gray shaded area is overplotted 
to indicate the range of ratios based on widely used synthesis models for solar metallicity 
and a Salpeter IMF with mass limits of 0.1 and 100 M yr _1 . These have been computed by 
Iglesias-Paramo et al. (2004) and Meurer et al. (2009) for the synthesis codes of Leitherer 
et al. (1999; Starbust99), Bruzual & Chariot (2003; BC03) and Fioc & Rocca-Volmerange 
(1997; PEGASE). All of the models primarily adopt stellar evolutionary tracks from the 
Padova group (e.g., Girardi et al. 1996 and references therein), but differ in their treatment 
of stellar atmospheres. There are uncertainties of ~20% in the Ha-to-FUV ratio due to the 
models alone. The dust-corrected ratios computed in the last section for the more luminous 
galaxies in our sample are within the range of expected model values. It thus appears 
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reasonable to use galaxies with SFR(Hct) > 3 x 10~ 2 M Q yr -1 to empirically define the 
fiducial Ha-to-FUV ratio at solar metallicity. Hereafter, we measure deviations from there 
instead of the K98 value. 



4.3. Metallicity 

Standard SFR conversion recipes (and hence the expected Ha-to-FUV ratio) generally 
assume solar metallicity populations. However, our sample spans a range of metallicities and 
is increasingly dominated by metal-poor dwarfs at low luminosities and SFRs. Given this 
trend in sample properties, it is possible that the variation in the Ha-to-FUV ratio could 
potentially be driven by systematic variations in metallicity. 

Metallicity influences the SED (spectral energy distribution) through its effect on the 
stellar opacity. Lower metallicity stars of a given mass will have lower opacities, lower 
pressures, and thus will be relatively smaller and have hotter atmospheres. They produce 
a larger number of UV photons (both ionizing and non-ionizing), so SFR conversions based 
on solar metallicity populations will tend to over-estimate the true SFR when applied to 
metal-poor systems (e.g., Lee et al. 2002, Brinchmann et al. 2004, Lee et al. 2009a). In the 
same vein however, metal-poor populations will produce more ionizing flux relative to non- 
ionizing UV continuum, leading to larger Ha-to-FUV ratios at low metallicity. This produces 
an effect which is the opposite of that observed, as also previously noted by Sullivan et al. 
(2000), Bell & Kennicutt (2001) and Meurer et al. (2009), and therefore cannot be the cause 
of the discrepancy. 

We illustrate these effects in Figure EJ which presents the SFR conversion factors and ex- 
pected Ha-to-FUV ratio, based on the Bruzual & Chariot (2003) stellar population synthesis 
models with Salpeter IMF and "Padova 1994" tracks. Values ranging from Z /5O (charac- 
teristic of the nebular oxygen abundance of the most metal poor local galaxies known; Kunth 
& Ostlin 2000) to 2.5Z are shown. Metallicity has less of an effect on the FUV than the Ha 
(note that the vertical scales on all panels span the same range), resulting in larger Ha-to- 
FUV ratios at low metallicity. Based on these models, the difference in the expected ratios at 
sub-solar metallicities relative to solar is small, <0.07 dex, so this should minimally impact 
the observed trends in Figure [51 We check on the magnitude of the effect by first estimating 
metallicities from the luminosity-metallicity relationship found locally for dwarf irregular 
galaxies (Lee et al. 2003), and then computing corrections for the SFRs of the individual 
galaxies in the sample. Replotting the metallicity-corrected SFRs confirms that the effect 
on the magnitude of the overall systematic decrease in the Ha-to-FUV ratio is negligible. It 
should be kept in mind, however, that the evolution of massive stars at sub-solar metallicity 
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is currently ill-constrained due to uncertainties in stellar rotation and mixing which may 
amplify this effect (e.g., Gallart et al. 2005, Leitherer 2008 and references therein). 

4.4. Ionizing Photon Loss 

In the calibration used to convert the Ha luminosity to the SFR, it is assumed that every 
Lyman continuum photon emitted results in the ionization of a hydrogen atom. However, 
if there is leakage of Lyman continuum photons into the IGM, the observed flux would be 
lower than expected, and the true SFR would be underestimated from Ha. In order for 
leakage alone to explain the observed trend in the Ha-to-FUV ratio, the fraction of ionizing 
photons escaping the galaxies would have to increase with decreasing luminosity and SFR, 
and at least half of the ionizing photons would have to be lost in the typical system with log 
L(Ha)~38.4. While the gas column densities in dwarf galaxies tend to be lower (e.g., Bigel 
et al. 2008), and there could be more leakage from the individual HII regions (the fraction 
of diffuse ionized gas appears to be nominally higher in dwarfs; e.g., Oey et al. 2007), star 
forming dwarf galaxies are often embedded in large envelopes of HI (e.g., Begum et al. 2008) 
which makes it unlikely that the Lyman continuum photons find their way completely out of 
a galaxy and into the IGM. Moreover, observations which have attempted to directly detect 
escaping Lyman continuum photons have all found upper-limits <10% (e.g., Leitherer et 
al. 1995, Bergvall et al. 2005, Grimes et al. 2007, Siana et al. 2007). These detection 
experiments have been performed on starbursting galaxies where leakage is thought to most 
likely occur. In contrast, the majority of the dwarfs in the sample are relatively normal in 
their current star formation activity (Lee et al. 2009a). It thus seems unlikely that leakage 
can be the cause for the trend. 

Ionizing photons can also be "lost" in a second way, via absorption of Lyman continuum 
photons by dust (i.e., dust absorption that Balmer decrement based attenution corrections 
cannot account for since the ionizing photons encounter dust before they are able to ionize 
hydrogen). However, it seems even more unlikely that this effect can cause the observed trend 
in the Ha-to-UV ratio since the absorption would need to be larger in relatively transparent 
dwarf galaxies than more luminous, dustier systems^ 



9 Further discussion on the direct absorption of Lyman continuum photons by dust in more luminous 
systems can be found in Buat et al. (2002), Hirashita et a. (2003) and Iglesias-Paramo et al. (2004). Such 
absorption is unlikely to be a salient issue in understanding the systematic decline in the Ha-to-UV ratio so 
we do not address it in any detail here. We do note however that the fraction of ionizing photons absorbed 
by dust cannot be very large, and is probably less than ^10-15%, based on the relative agreement of dust 
corrected UV and Ha SFRs at the luminous end of our sample. 
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4.5. Starbursts in Dwarf Galaxies 

Previous work tentatively attributed possible trends in the Ha-to-FUV ratio to a sys- 
tematic increase in the prevalence of bursts in the recent star formation histories of low mass 
galaxies (e.g., Bell & Kennicutt 2001). The fiducial expected ratio is an "equilibrium" value 
which is calculated assuming that the SFR has been constant over a time period that is 
long enough for the birth of stars responsible for the FUV and Ha emission to balance their 
deaths. Variations in the SFR over timescales ~100 Myr would disrupt this equilibrium. In 
the time following a burst of star formation, a deficiency of ionizing stars develops as they 
expire relative to lower-mass, longer-lived stars that also contribute to the UV emission, 
and the Ha-to-FUV flux ratio will thus be lower than expected (e.g., Sullivan et al. 2000, 
Sullivan et al. 2004, Iglesias-Paramo et al. 2004). 

To examine whether variations in the SFR can account for the decline of the Ha-to- 
FUV ratio of the magnitude that is observed, synthesis modeling is needed. Here the work 
of Iglesias-Paramo et al. (2004) can be directly applied. Iglesias-Paramo et al. (2004) sought 
to use the Ha-to-FUV flux ratio to constrain the effects of the cluster environment on the 
recent SFH of spiral galaxies, and computed a grid of Starburst99 models for a range of burst 
amplitudes and durations, assuming a Salpeter IMF. The bursts were superimposed on 13 
Gyr-old models which were evolved assuming an exponentially declining SFR, with decay 
timescales typical of values normally used to fit star forming galaxies (from 3 to 15 Gyr). 
Variations in the predicted Ha-to-FUV flux ratio were found to be insensitive to differences 
in decay timescale values within this range. The results of their modeling are tabulated, 
and it is shown that in order to produce Ha-to-FUV flux ratios that are on average about a 
factor of two lower than the fiducial equilibrium value in a given sample, bursts which elevate 
the SFR of the galaxies by a factor of 100 for a 100 Myr duration are required. However, 
such large amplitude bursts do not appear to commonly occur in the overall dwarf galaxy 
population. The most direct constraints on burst amplitudes and durations in dwarf galaxies 
are provided by studies which reconstruct star formation histories from resolved observations 
of stellar populations in the nearest low-mass systems. Most recently, Weisz et al. (2008) and 
McQuinn et al. (2009) have analyzed HST/ACS imaging of dwarf galaxies in the M81 group 
and a small sample of low-mass starburst systems, respectively. Typical burst amplitudes 
range from a few to ~10, an order of magnitude smaller than the factor of 100 bursts required 
to depress the Ha-to-FUV flux ratios by a factor of 2. Further, the 11HUGS sample itself 
provides a statistical constraint on the average dwarf galaxy starburst amplitude via the 
ratio of fraction of star formation (as traced by Ha) concentrated in starbursts to their 
number fraction (Lee et al. 2009a). The high degree of statistical completeness of 11HUGS 
makes this calculation possible for galaxies with Mb < —15, and again, the burst amplitude 
is found to be relatively modest (~4), at least for this population of dwarfs. From the above 
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arguments it appears that modest starbursts by themselves, even if they are common in 
low-mass systems, cannot explain the observed deviation in the Ha-to-FUV flux ratios from 
the expected value. 

4.6. Stochasticity in High Mass Star Formation at Low SFRs 

An implicit assumption in all stellar population synthesis models (and the resulting 
calibrations based upon them) is that statistical numbers of stars spanning a full range of 
masses, typically to a limit of M*= 100 M Q , are produced. In the regime of ultra- low SFRs 
where only a handful of O-stars are formed in a given system over timescales comparable 
to their lifetimes (i.e., a few million years), this assumption is clearly violated. Under these 
circumstances, Ha emission can appear deficient, or even absent, although lower mass star 
formation, which dominates the total mass formed, may be taking place and the underlying 
IMF is invariant. Such statistical effects on the Ha luminosity has been previously explored 
for example by Boissier et al. (2007) and Thilker et al. (2007) in the context of extended 
UV disks. More detailed work has been done by Cervino & Valls-Gabaud (2003), Cervino 
et al. (2003) and Cervino & Luridiana (2004) for modeling properties of low mass stellar 
clusters and other simple, single-age stellar populations. To evaluate whether such effects 
can account for our observations, we first perform a simple back-of-the-envelope calculation, 
and then examine results from the Monte Carlo simulations of Tremonti et al. (2009). 

We first estimate when Poisson fluctuations should begin to affect the Ha output by 
computing the SFR at which the number of O-stars on the main sequence at any given time 
will be less than 10. Assuming a Salpeter IMF with standard mass limits of 0.1-100M o yr _1 , 
and considering that only stars with masses >18 M significantly contribute to the ionizing 
flux (e.g., Hunter & Elmegreen 2004, appendix B), simple integrals show that there are 10 
ionizing stars for a total stellar mass of 4.3xl0 3 M Q formed. To compute the corresponding 
SFR above which the Ha flux should be robust against stochasticity in the formation of 
ionizing stars, we divide this mass by the 3 Myr lifetime of a 100 M Q star (Schaerer et al. 
1993), which yields 1.4xl0 -3 M Q yr _1 (log SFR=-2.8). This is a conservative limit compared 
to calculations based on the formation of single O-stars (e.g., Lee et al. 2009a; Meurer et al. 
20090 However, Figures [1] and [5] show that systematic deficiencies in the Ha-to-FUV ratio 
begin at SFRs that are almost an order of magnitude higher (log SFR~-2.0). At log SFR=- 



10 In Lee et al. (2009a), the calculation at the end of §2.1.1 was intended to estimate the SFR required to 
produce a single O-star at any given time for a Salpeter IMF with mass limits of 0.1 and 100 Mq, and not 
10 O-stars as stated there. 



-19- 



2.8 the ratio is not just beginning to deviate from the expected value; it is already more 
than a factor of two lower. This back-of-the envelope calculation suggests that stochasticity 
in the sampling of an invariant Salpeter IMF alone is unlikely to give rise to the observed 
trends. The same holds true for the Kroupa (2001) and Chabrier IMFs (2003), for which 
the computed SFR limit would be 1.5 times lower. In all cases the FUV flux is negligibly 
affected. 

Tremonti et al. (2009) have preformed detailed Monte Carlo simulations to confirm 
these conclusions. Such calculations not only allow for a more robust prediction of the mean 
deviations due to stochasticity in the production of massive ionizing stars, but also provide 
quantitative predictions of the accompanying increase in scatter. In brief, we generate a 
pool of stars which have a mass distribution given by the Salpeter IMF, and assign each 
star a mass- dependent main sequence lifetime, ionizing photon production rate, and FUV 
luminosity. We randomly draw stars from the pool at a rate proportional to the desired 
SFR and follow the ensemble population over time. SFRs are calculated from the FUV and 
Ha luminosities at t>l Gyr. We have run simulations for a large grid of SFRs; the results 
are overplotted on the data in Figure CD The median predicted values (solid line) as well as 
the values at the 2.5 and 97.5 percentile points (dotted lines) are shown. The median flux 
ratio begins to significantly deviate from the fiducial value only for log SFR< —3, which is 
consistent with expectations based on the back-of-the-envelope estimate given above. 

We note that the results of our calculations agree with the work of Cervino et al. (2003) 
and Cervino & Luridiana (2004), who have used statistical formalism to compute values for 
the "Lowest Luminosity Limit" (LLL), the limiting stellar mass formed in a single event 
whose properties can be modeled without accounting for stochastic effects. According to 
their calculations, ~3xl0 3 M is the mass associated with the LLL in the production of 
Lyman continuum photons. Again, taking the O-star lifetime to be ~3 Myr as above, this 
limiting mass would correspond to log SFR=-3. However, as Cervino et al. (2003) discuss, 
statistical effects may be present in formation events involving up to ten times the LLL 
mass. Based upon the Monte Carlo simulations, these effects manifest as an increase in 
scatter in the Ha-to-FUV ratio beginning at log SFR~-2, although the median value of the 
ratio remains stable until log SFR~-3. 

Therefore, to summarize, although stochasticity in the sampling of an invariant IMF 
does result in a declining Ha-to-FUV ratio with SFR, we find that the mean effect is not 
large enough to account for the observations. One important caveat regarding this conclusion 
however is that our calculations assume that the SFR is constant. As discussed in Section I4T5] 
moderate variations of the SFR with time are seen in the recent star formation histories 
reconstructed from resolved stellar population data of nearby dwarf galaxies (e.g., Dohm- 
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Palmer et al. 1998; Weisz et al. 2008), and stochasticity may amplify the effect of such 
variations on the Ha-to-FUV ratio. As illustrated in Figure [7] however, an initial analysis 
where bursts with 100 Myr duration and factor of four amplitude are included in our Monte 
Carlo simulations show that this is unlikely to produce significant changes in the predictions 
based on a constant SFR alone. However, more comprehensive modeling covering a range of 
star formation histories is needed, and these issues are further examined in Tremonti et al. 
(2009). In particular, histories where the SFR shows short timescale (~10 Myr) fluctuations 
of factors of 10 or greater may be considered. Such histories may be characteristic of the 
most extreme dwarfs in our sample, in which only one or two HII regions are visible, and 
where such fluctuations are not necessarily indicative of starburst activity, but rather of 
Poisson noise in the formation of individual star clusters. 



4.7. IMF 

Finally, we examine the basal assumption required for the calibration of all star for- 
mation indicators, the form of the IMF. Naturally, indicators which trace highly luminous 
but relatively rarer massive stars are particularly sensitive to its form. However, the IMF is 
generally assumed to be invariant, and if so, the Ha-to-FUV flux ratio should be constant, 
modulo the other factors discussed above. Uncertainties in a universal IMF therefore cannot 
themselves produce trends in the Ha-to-FUV flux ratio. Rather, systematic variations which 
result in a deficiency of ionizing stars in dwarf galaxies would be needed to explain the ob- 
servations. This possibility must be considered, given that we cannot explain the magnitude 
of the Ha-to-FUV offset with the other parameters that have been explored thus far. 

Evaluating the plausibility of this scenario in present day galaxies has been an issue 
of much recent work. For example, Hoversten & Glazebrook (2008; HG08) examined the 
Ha EW (equivalent width) and optical colors of a large sample of galaxies taken from the 
Sloan Digital Sky Survey (SDSS). They find that the Ha EWs are systematically low for 
the optical colors in low-luminosity galaxies given expectations based on continuous star 
formation. This is consistent with our results for a more limited sample drawn from the 
Local Volume (Lee 2006, 2008). However, whereas Lee (2006, 2008) concluded that such 
discrepancies could be explained by invoking starbursts in the recent past histories of dwarfs 
(most systems would then be in a post-burst state where the ionizing stars from the burst 
have died off, but the intermediate mass stars continue to contribute to the continuum flux 
density at Ha, thus depressing the Ha EW), the analysis of HG08 appears to have ruled 
this out. The much larger SDSS sample of HG08 enabled more detailed modeling which 
probed the temporal variation of the EWs and colors over a burst cycle. HG08 found that 



21 



fine-tuning to coordinate burst times was necessary to reproduce the data, and that a more 
likely explanation was that low-luminosity galaxies have an IMF deficient in massive stars. 

Recently, Meurer et al. (2009) have examined correlations in the Ha-to-FUV flux ratio 
with the Ha and i?-band surface brightnesses, using the Hl-selected samples of the Survey 
of Ionization in Neutral Gas Galaxies (SINGG; Meurer et al. 2006) and the Survey of UV 
emission in Neutral Gas Galaxies (SUNGG). Meurer et al. find that the Ha-to-FUV flux 
ratio decreases with decreasing surface brightness. After ruling out the more likely causes 
of the systematic in an analysis parallel to the one presented here, they have concluded 
that it is evidence for a non-uniform IMF. Their study is complementary to ours since they 
have examined correlations in the Ha-to-FUV flux ratio with the Ha and i?-band surface 
brightness instead of integrated luminosities, and use an independent dataset. The galaxy 
sample however is weighted more heavily toward more luminous star-forming systems (e.g. 
Paper I, Figure 7 & 8); there are only ~10 galaxies with dust corrected SFR(Ha)<0.01 
M Q yr _1 and this precludes a robust determination of the observational trend with integrated 
luminosity as presented here (M. Seibert, private communication). 

On the theoretical side, Kroupa & Weidner (2003) and Weidner & Kroupa (2005, 2006) 
have formulated a model based on both statistical arguments and observational constraints 
which produces IMFs that appear to be steeper for galaxies with lower SFRs. Several critical 
assumptions are made. The first is that all stars are born in stellar clusters, which universally 
form stars according to the Kroupa (2001) IMF. The formation of stellar clusters themselves 
is also governed by a power law. While the convolution of cluster and stellar IMFs by itself 
does not produce effects beyond simple random sampling of the stellar IMF alone, the model 
also assumes that (i) the maximum mass of a cluster formed in a given galaxy, m™ ax , depends 
on its SFR, and (it) the maximum mass of a star formed in a given cluster, m™ ax depends 
on the total cluster mass M d , and is lower than the strict limit where m™ ax = M d . The 
claim is that both m d ax (SFR) and m™ ax (M d ) are both robustly constrained by observations. 
Further, it is shown that m™ ax (M c i) is also well-modeled by the statistical condition: 





All these effects taken together can be expressed as: 




(9) 



where ^igimf is referred to as the "Integrated Galaxial IMF. 
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Recently, Pflamm-Altenburg, Weidner & Kroupa (2007; 2009) have used this model to 
generate predictions for trends in FUV and Ha properties. They consider the case where the 
cluster mass function has a Salpeter slope (referred to as the "standard model" ) and also 
examine the consequences of assuming one which has a higher ratio of high-mass to low-mass 
clusters, i.e., a d =1.0 for 5M < M d < 5OM and a d =2.0 for 50M Q < M cl < m^ ax (SFR). 
Adopting this £ c j results in smaller relative deviations from a universal Kroupa IMF and 
is referred to as the "minimall" model. The predictions are overplotted on the data in 
Figure [HI where the open red squares represent the standard model, while the filled red 
squares represent the minimall model. Both are in good agreement with the data at the 
faint end, while the data for the more luminous galaxies in the sample prefer the minimall 
model. Therefore, the IGIMF model is able to account for the declining Ha-to-FUV ratio 
observed in our Local Volume galaxy sample. 

A next critical test of the IGIMF model would involve computing the expected scatter 
in the relationship given plausible star formation histories, since it is not clear whether 
the IGIMF will allow for Ha-to-UV ratios higher than the predicted curve (as seen in the 
data, modulo extinction), given the deterministic upper- mass limits on the IMF in the low 
activity regime. Monte Carlo simulations analogous to the ones described in the previous 
section should also be run with the IGIMF as the input IMF to check that the combination 
of stochasticity, variations in the SFH and the IGIMF do not lead to larger declines in 
the Ha-to-UV ratio than is seen in the data. Finally, further work is still also required to 
understand the relationship between the IGIMF model and the physics of the fragmentation 
and collapse of gas into stars. 



4.8. L(Hct) as an SFR Indicator at Low Activities 

Whether or not IMF variations (or other unidentified issues) are the true culprits under- 
lying the systematic variation in the Ha-to-FUV ratio, it seems reasonable to surmise that 
the FUV luminosity is the more robust SFR indicator in individual galaxies with low total 
SFRs and low dust attenuations. This is also supported by the finding that the increase 
in scatter at low galaxy luminosities in the Ha equivalent width (which traces the SFR per 
unit stellar mass), as discussed in Lee et al. (2007), is mitigated when SFR per unit mass 
is instead calculated from the FUV. Figure M illustrates the variation of the specific SFR 
(SSFR; the SFR per unit stellar mass) with M B where the SSFR has been dust corrected and 
is based on the Ha (left panel) and FUV luminosities (right panel). The SFRs and attenua- 
tions corrections are identical to those used throughout this paper, while the stellar masses 
are based upon Mb and color-dependent mass-to-light ratios as described in Bothwell et al. 
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(2009). These plots are analogous to those presented in Lee et al. (2007). The previously 
reported increase in scatter near Mb ~ —15 is apparent in both panels. However, whereas 
the scatter in SSFR(Ha) increases toward the lowest luminosities by ~60%, the increase is 
only ~25% for SSFR(FUV). Given the longer lifetimes and greater numbers of stars which 
dominate the UV output of a galaxy, this result is reasonable as the UV will be less affected 
by purely stochastic variations in the formation of high mass stars. This relative reduction 
in the increase in scatter is also generally consistent with the predictions of our Monte Carlo 
simulations for SFR<0.01 M Q yr _1 as reported in Tremonti et al. (2009) and discussed in 
Section SH 

Therefore, the results of our analysis suggest that when possible, the SFR should be 
based on the FUV luminosity when SFR<0.01 M yr -1 . Ha is a less reliable indicator when 
the integrated, galaxy-wide log L(Ha) [ergs s _1 /< 39.4. However, for those instances where 
only an Ha measurement is available, we provide an empirical re-calibration of SFR(Ha), 
using the extinction corrected K98 SFR(FUV) as the reference value. For an observed, 
non-dust corrected, integrated L(Ha)< 2.5 x 10 39 ergs s" 1 : 

log (SFR [Meyr^ 1 ]) = 0.62 log(7.9 x 10" 42 L(iJa)[ergs s" 1 ]) - 0.47. (10) 

This relation is the result of a linear least squares fit of SFR(FUV), where L(FUV) has 
been corrected for internal dust extinction, as a function of SFR(Ha), which has not been 
dust corrected. Thus, the SFR calculated using this expression will represent intrinsic, 
dust unattenuated values as appropriate for the dust content typical of local galaxies, trace 
the activity averaged over the timescale of the FUV emission (~100 Myrs) instead of the 
native Ha instantaneous timescale (~3 Myrs), and carries the same assumptions (fiducial 
stellar models, Salpeter IMF) as the K98 UV SFR calibration. The random error on the 
SFR, which will increase with decreasing L(Ha), can be estimated from the l-a scatter in 
SFR(Ha)/SFR(FUV) as listed in Table 2. 

As with any SFR prescription, it is important to be aware of the ranges of applicability 
of Equation [TD] and its limitations. The relation has been calibrated with integrated SFRs of 
local dwarf irregular galaxies for which log L(Ha)>37, so application to other systems with 
low total SFRs (e.g., massive elliptical galaxies, or galaxies with log L(Ha) [ergs s _1 ]<37) 
may not be warranted!"] We also note that higher order terms may be required to more 
accurately describe the systematic underestimation of the SFR in this regime, although the 



11 Note that although the dataset extends to log L(Ha) [ergs s _1 ]~36, there are only 5 galaxies with log 
L(Ha) [ergs s _1 ]<37. Clearly, a larger sample is needed to constrain systematics in the Ha-to-FUV ratio 
and the empirical re-calibration in this regime. 
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current dataset does not adequately constrain such terms. For example, Pflamm-Altenburg 
et al. (2007) offer a calibration based on the IGIMF model that is described by a 5th order 
polynomial. In Figure [TUl our empirical relation is compared with the IGIMF prescriptions. 
The IGIMF "minimall" model is consistent for log L(Ha)>38, but predicts up to 0.2 dex 
more star formation in the regime where there is sufficient data available to calibrate our 
relation. 



5. Summary and Conclusions 

Using Ha flux measurements and GALEX UV photometry from the 11HUGS Legacy 
program (Kennicutt et al. 2008; Lee et al. 2009b), we evaluate the consistency between star 
formation rates (SFRs) inferred from Ha and the far UV (FUV) continuum over 5 orders 
of magnitude, down to ultra-low SFRs of ~0.0001 M Q yr _1 , assuming standard conversion 
recipes in which the SFR scales linearly with luminosity at a given wavelength. The primary 
observational result is that the Ha luminosity underestimates the SFR relative to the FUV 
luminosity in dwarf galaxies that are roughly less active than the Small Magellanic Cloud; 
that is, the observed Ha-to-FUV flux ratio systematically decreases with declining lumi- 
nosity. The effects of uncertainties in the stellar evolution tracks and model atmospheres, 
non-solar metallicities, non-constant SFHs, possible leakage of ionizing photons/departures 
from Case B recombination, dust attenuation, stochasticity in the formation of high-mass 
stars, are all considered. However, none of these potential drivers acting alone are able to 
explain the magnitude of the observed systematic. The underlying cause for the trend is not 
clear, although we cannot rule out variations in the IMF which result in a deficiency of high 
mass, ionizing stars in dwarf galaxies. In fact, the semi-empirical IGIMF model of Kroupa 
Sz Weidner (2003), which results in such a deficiency, predicts a trend in the Ha-to-FUV flux 
ratio that is strikingly consistent with our data. It is also possible that some combination of 
effects may conspire to produce the observed systematic. 

More work is needed to further elucidate the nature of the systematic as presented in 
this paper and in Meurer et al. (2009). For example, stochasticity may amplify the effects 
of bursty or non-uniform "flickering" star formation histories on the Ha-to-FUV ratio, as 
discussed in Section 14.61 and we investigate such effects in greater detail in Tremonti et 
al. (2009). Another obvious issue is the difficulty in evaluating whether the perceived 
inconsistency is simply due to unknowns in stellar evolution modeling of massive stars at 
low metallicity. Also, although it seems unlikely that the leakage of ionizing photons can 
be significant enough to account for the low Ha-to-FUV flux ratio observed, it remains a 
possibility that photons can preferentially escape in a direction perpendicular to the disk 
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in low mass galaxies. Tracing the formation of high mass, ionizing stars in an independent 
way, such as with radio continuum measurements, may be useful. Examining the relative 
spatial distributions of emission within UV-bright, Ha-faint galaxies may also provide clues 
on the origin of the systematic. In a number of dwarfs, we observe that whereas the FUV 
emission extends over the optical disk, the Ha (and infrared emission) is far more localized 
in a few scattered clumps. Finally, to further test the possibility of systematic variations in 
the IMF, other consequences must be followed up to ensure that at a minimum, these yield 
a consistent picture. As pointed out by Koppen et al. (2007), IMFs deficient in high mass 
stars should imprint signatures in the observed metal abundance ratios, and these must be 
carefully reexamined in this context. 
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Table I. 11HUGS Ha and UV SFRs 



# 


Galaxy name 


RA 


DEC 


E(B - V) 


D 


method 


M B 


T 


A(Ha) 


A(FUV) 


code 


SFR(UV) 


SFR(Ha)/SFR(UV) 






[J2000] 


[J2000] 


[mag] 


[Mpc] 




[mag] 




[mag] 


[mag] 




LOG [Mq yr- 1 ] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


1 


WLM 


000158.1 


-152739 


0.04 


0.92 


trgb 


-13.94 


10 


0.10 






-2.21 


0.31 


2 


ESO409-IG015 


000531.8 


-280553 


0.02 


10.4 


v(fiow) 


-15.0 


6 


0.14 






-1.45 


1.48 


3 


NGC24 


000956.7 


-245744 


0.02 


8.1 


v(fiow) 


-17.4 


5 


0.41 


0.74 




-0.70 


0.56 


4 


NGC45 


001404.0 


-231055 


0.02 


7.1 


v(fiow) 


-18.01 


8 


0.50 


0.30 




-0.40 


0.99 


5 


NGC55 


001454.0 


-391149 


0.01 


2.17 


trgb 


-18.32 


9 


0.55 


0.86 




-0.26 


0.86 



Note. — Col.(l): Running index number in this table. Col.(2): Galaxy name. Col. (3-4): J2000 R.A. and Dec. from NED. Col. (5): Milky Way reddening based on 
the maps of Schlegel et al. (1998). Col. (6): Distance to the galaxy in Mpc. Col. (7): Method of distance determination as compiled, described, and referenced in Paper 
I. H =75 km s _1 Mpc -1 is assumed for distances estimated from Local Group flow corrected recessional velocities. Col. (8): Absolute B-band magnitude based on the 
photometry compiled in Paper I, and the distances and Milky Way reddenings given the previous columns. Col. (9): RC3 Morphological T-type, compiled as described 
in Paper I. Col. (10): Ha nebular attentuation in magnitudes calculated as described in Section . Col. (11): FUV attentuation in magnitudes based on the TIR/FUV 
ratio as described in Section . Col. (12): Source of the H«/H/3 Balmer line ratio used to calculate the Ha nebular attentuation. "mk" refers to the integrated spectral 
atlas of Moustakas and Kcnnicutt (2006), and "nfgs" refers to the integrated spectra of the Nearby Field Galaxy Survey (Jansen et al. 2000). A dash preceding the 
reference ("-mk" or "-nfgs") indicates that Hct/H/3 is less than the expected Case B ratio of 2.86. This results in a negative extinction correction, so in these cases the 
value of A(Ho) given is estimated from Mb, as are those where spectral line ratio measurements arc not available from one of these sources. Col. (13): SFR, in units of 
log[M0 yr -1 ], based on the conversion of K98, the integrated GALEX FUV magnitude reported in Paper II and the A(FUV) given in this table. If no value for A(FUV) 
is listed, the attenuation is estimated from 1.8*A(Ha). Col. (14): Ratio of SFRs calculated from H« and the FUV. The Ho SFRs arc based on the conversion of K98, 
line fluxes and [NIIJHct corrections given in Paper I, and the A(Ha) listed in this table. 
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Table 2. Ha-to-UV SFR Ratios 



Observed 



LOG SFR(Ha) 


j SFR(Hq) 
SFR(UV) 


lcr-scatter 


Ngal 


0.25 


0.20 


0.37 


11 


-0.25 


0.17 


0.30 


21 


-0.75 


0.07 


0.26 


39 


1 

- 1 ,zo 


n no 


U.zo 


DO 


-1.75 


-0.10 


0.22 


55 


-2.25 


-0.23 


0.22 


61 


-2.75 


-0.46 


0.26 


33 


-3.50 


-0.49 


0.58 


20 


-4.50 


-1.29 


0.57 


6 


M B 


Trir , SFR(Hq) 
SFR(UV) 


l<T-scatter 


Ngal 


-20.5 


0.18 


0.23 


14 


-19.5 


0.20 


0.39 


12 


-18.5 


0.08 


0.34 


27 


-17.5 


-0.01 


0.31 


36 


-16.5 


0.00 


0.21 


35 


-15.5 


-0.22 


0.37 


58 


-14.5 


-0.26 


0.29 


61 


-13.5 


-0.31 


0.30 


41 


-12.5 


-0.41 


0.66 


21 


-11.5 


-0.48 


0.61 


7 


Dust Corrected 


LOG SFR(Ha) 


T ^ SFR(Hq) 
LOG SFR(UV) 


la-scatter 


Ngal 


0.50 


-0.14 


0.29 


20 


-0.25 


-0.09 


0.20 


26 


-0.75 


-0.13 


0.22 


44 


-1.25 


-0.13 


0.22 


51 


-1.75 


-0.18 


0.17 


58 


-2.25 


-0.27 


0.22 


56 


-2.75 


-0.51 


0.25 


31 


-3.50 


-0.55 


0.57 


21 


-4.50 


-1.43 


0.59 


5 


M B 


yFR(Ha) 

SFR(UV) 


lcr-scatter 


N gal 


-20.5 


-0.21 


0.23 


14 
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-0.10 


0.36 


12 


-18.5 
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0.23 


27 


-17.5 


-0.16 


0.17 


36 


-16.5 


-0.12 


0.18 


35 


-15.5 


-0.27 


0.36 


58 


-14.5 


-0.30 


0.28 


61 


-13.5 


-0.34 


0.30 


41 


-12.5 


-0.45 


0.67 


21 


-11.5 


-0.51 


0.61 


7 
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LOG l_(Ha) [ergs s" 1 ] 




LOG SFR(Ha) [M yr" 1 ] 



Fig. 1. — Comparison of observed FUV and Ha SFRs (i.e., without corrections for internal 
dust attenuation), calculated using the linear conversion recipes given in Kennicutt (1998). 
Axes that indicate the corresponding Ha and FUV luminosities are also shown. Morpho- 
logical type is distinguished with different symbols as indicated. The solid line represents 
a one-to-one correspondence between the SFRs, while the dashed line is an ordinary least 
squares bisector fit to the data, given by log[SFR(FUV)]=0.79 log[SFR(Ha)] - 0.20. The 
black filled symbols represent galaxies where the Ha flux may be underestimated because 
the narrowband imaging did not wholly enclose the galaxy and are excluded from the fit. 



-33- 



LOG [L(Ha) ergs s" 1 ] 
36 38 40 



42 



1 



1 

2 



T 



TT 



b t 't .. fiti&gm j 



t — i — i — i — i — =n 

NGC3593 NGC3034 - 



:. ® . . -r- 



IDDO210 



1 - — 



UGC2847 



UGCA4' 

S 

J I I L 



1 



-4 



-2 







LOG [SFR(Ha) M yr" 1 ] 



-5— t 



M 1 1 1 I 1 1 I 

KDG61 

1 E- 



o 
l 

2 



I 1 1 1 I 1 

NGC3034 



NGC3593 



E_L 



. UGC8215 
DD02 iG6A438 
I 



J i I i L 



UGC2847-I 



J I I I = 



-10 -12 -14 -16 -18 -20 -22 



Fig. 2. — The ratio of Ha-to-FUV SFRs, where corrections for internal dust attenuation have 
not yet been applied (or equivalently, the observed Ha-to-FUV flux ratio), as a function of 
the Ha SFR/luminosity (top panel) and the B-band absolute magnitude (bottom panel). 
The solid lines indicate ratios of unity, while the dashed line shows the linear least squares 
fits to the data, which are given by log[SFR(Ha)/SFR(FUV)]=0.26 log[SFR(Ha)]+0.30, and 
log[SFR(Ha)/SFR(FUV)]=-0.08 M B -1.39. Binned averages of the data with la scatter are 
also overplotted (blue symbols) and are listed in Table 2. Circled points represent galaxies 
where the Ha flux may be underestimated because the narrowband imaging did not wholly 
enclose the galaxy, and are excluded from the fit and statistics. 
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Fig. 3. — A(Ha), the Ha nebular dust attenuation in magnitudes, based on measurements 
of the Balmer decrement, plotted against the B-band absolute magnitude Mb for galaxies in 
the integrated spectral atlas of Moustakas & Kennicutt (2006). The black filled circles are 
star-forming galaxies, and those with A(Ha) > have been included in the over-plotted fit 
(see Equation [5]) . The various crosses represent galaxies which have been excluded from the 
fit because they have line ratios indicating the presence of an AGN (red), signal-to-noise in 
H/? <10 (blue), or ¥L(3 equivalent width < 5A (magenta). Luminous star-forming galaxies 
exhibit greater levels of nebular attenuation than those at lower luminosity. The over-plotted 
fit is used to estimate the nebular extinction when measurements of the Balmer decrement 
from integrated spectral measurements are not available. 



-35 - 




Fig. 4. — A comparison of different estimates of the Ha and FUV attenuation in galaxies 
in the If HUGS sample. Bottom panel: FUV attenuations, A(FUV), based on T1R/FUV 
ratios computed from MIPS photometry, and Ha attenuations A(Ha) based on (i) integrated 
spectroscopic measurements of the Balmer decrement (solid large points) and, (ii) A(Ha) 
estimated from a Mb dependent scaling relationship calibrated with Balmer decrement mea- 
surements (open points). Top panel: A(FUV) based on the linear IRX-/3 relationship of 
Cortese et al. (2006) and our best available value for A(Ha) (small black points). These 
latter A(FUV) estimates appear to be too high and are not used in our analysis. The solid 
line in both panels represents A(FUV) = f .8A(Ha), as expected for the Calzetti obscuration 
curve and differential extinction law. The Calzetti law provides a good description of our 
most robust attenuation measurements (solid large points). 
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Fig. 5. — Same as Figure El but corrected for attenuation by dust internal to the galax- 
ies as described in Section [TJ Data points which have nebular attenuations computed 
from the Balmer decrement are distinguished by red X's, and those where the FUV at- 
tenuations are based on the TIR-to-FUV ratio are shown as red crosses. Galaxies for 
which both Balmer decrement and TIR-to-FUV based corrections are applied thus ap- 
pear as red stars. Corrections for the remaining points are estimated with scaling rela- 
tionships. The shaded band represents the range of Ha-to-FUV ratios predicted by com- 
monly used stellar population models as described in Section 14.21 The dashed lines in 
the top panel are given by log[SFR(Ha)/SFR(FUV)]=-0.13 for log[SFR(Ha)]>-1.5, and 
log[SFR(Ha)/SFR(FUV)]=0.32 log[SFR(Ha)]+0.37 for log[SFR(H«)]<-1.5. In the bottom 
panel the fit is given by log[SFR(Ha)/SFR(FUV)]=-0.05 M B -0.99. 
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Fig. 6. — Star formation rate conversion factors and the expected Ha-to-FUV ratio as a 
function of metallicity calculated from the Bruzual & Chariot (2003) population synthesis 
models assuming a Salpeter IMF. The top panel shows the Ha conversion factor in units 
of log (M yr _1 )/(ergs s _1 ), assuming case B recombination with nebular temperatures and 
densities of 10 4 K and 100 cm -3 . The middle panel shows the conversion factor for the FUV 
flux density in units of log(M yr _1 ) / (ergs s _1 Hz), calculated using the GALEX FUV filter 
profile. The bottom panel shows the corresponding expected ratios. The solid horizontal line 
in each panel indicates the values as given in K98. In the bottom panel, the shaded region 
shows the range of ratios based on widely-used synthesis models as described in Section 14.21 
and also plotted in Figure HI The Ha-to-FUV flux ratio increases with decreasing metallicity. 
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Fig. 7. — Top panel: Same as Figure [U but corrected for attenuation by dust internal to 
the galaxies as described in Section [H and with predictions from the stochastic models of 
Tremonti et al. (2007; 2009) overplotted. Bottom panel: Same as the top panel of Figure El 
but again with the stochastic model predictions overplotted. The median predicted values 
(solid line) as well as the values at the 2.5 and 97.5 percentile points (dotted lines) are shown. 
The green curves represent a model where the star formation rate is constant at all times, 
while the red curves show a model where bursts with an amplitude of four and a duration 
of 100 Myrs are added every 500 Myrs. 
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Fig. 8. — Top panel: Same as Figure [U but corrected for attenuation by dust internal to 
the galaxies as described in Section [TJ Predictions based on the IGIMF model (Kroupa 
& Weidner 2003; Pflamm-Altenburg, Weidner & Kroupa 2009) are overplotted. Bottom 
panel: Same as the top panel of Figure [5j but again with predictions based on the IGIMF 
model overplotted. The open and filled red squares represent the IGIMF "standard" and 
"minimall" models, respectively. The bottom and left hand axes are shown in gray to signify 
that the K98 SFR scales would not be valid at low SFRs under the assumptions of the IGIMF 
model. 
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Fig. 9. — The variation of the specific SFR (SSFR; the SFR per unit stellar mass) with 
M B where the SSFR has been dust corrected and is based on the Ha (left panel) and FUV 
luminosities (right panel). These plots are analogous to those presented in Lee et al. (2007) 
which directly used the Ha EW to trace the SSFR. In that paper, a transition at Mb ~ — 15, 
characterized by an increase in scatter toward low luminosities, was identified. The same 
structure is apparent in these figures, but with the notable distinction that the relative 
increase in scatter is reduced when the SSFR is calculated from the FUV luminosity. 
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Fig. 10. — The empirical calibration of L(Ha) as a SFR indicator in the regime of low inte- 
grated activities based on the 11HUGS dataset (equation [TDJ solid line), compared with the 
Pflamm-Altenburg, Weidner & Kroupa (2007) predictions for the IGIMF model. The dotted 
and dashed curves represent the IGIMF "standard" and "minimall" models, respectively. 



